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II  I 


ABSTRACT 


The  amplitude  and  period  ol  over  600  long-period  vertical-component 

seismic  signals  were  measured  at  regular  intervals  in  order  to  characterize 

the  LK  signal  and  its  coda.  Signal  envelopes  out  to  20  minutes  after  LK  onset 

were  found  to  vary  considerably,  but  coda  amplitudes  decayed  rather  uniformly 

with  time.  Using  average  signal  and  coda  parameters  and  randomly-generated 

events,  a program  which  simulates  long-period  signal  i rcerf erence  at  seismic 

stations  was  devised  to  predict  interference  effects  on  network  detection 

statistics.  Increasing  the  size  of  a network  should  not  increase  the  number 

of  reported  events  with  interference  at  any  M , but  lowering  the  station 

thresholds  will  increase  tills  number  at  lower  M values. 
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INTRODUCTION 


space-t  lme  distribution  of  seismic  „ctIv,ty  ovcr  the  8lobe  ^ 
rcsu  ts  In  the  simultaneous  arrival  of  two  or  more  signals  at  a particular 
station.  This  phenomenon  is  especially  serious  in  a detection  context  when 
one  event  is  large  and  effectively  masks  the  smaller  arriving  signals  of  on, 
or  more  other  events.  The  amount  of  actual  interference  on  long-period 
s;at;o„s  has  been  documented  by  Mack  and  Robertson  (1,73)  and  von  Seggern 
, can  affect  10  to  20  percent  of  the  seismic  events  which  occur, 
an  so  in  detection  it  is  a significant  hindrance  which  is  not  accounted’ 
for  in  routine  prediction  of  seismic  network  capability.  Kor  „ net„ork.  one 
might  presume  that  the  addition  of  more  stations  at  various  sites  on  the 
globe  should  ameliorate  the  Intel ference  problem  since  there  are  then  more 
C ances  or  separation  in  arrival  time  of  signals  from  events  which  occur 
C ose  n time  but  not  in  space.  On  the  other  hand  the  interference  problem 
»ay  increase  in  severity  as  the  threshold  of  network  drops,  due  to  more 

rr  HSl"Ce  “ny  "°re  eVC"tS  b“~  ‘—ivi.  »ne  might  expect  though 
tha  the  advantage  of  a lower  threshold  outweighs  any  increased  interference 

problem.,  These  same  presumptions  apply  when  the  network  threshold  is 

decreased  by  improvements  in  existing  sites;  that  is,  by  the  lowering  of 
noise  background  on  recordings. 

niese  presumptions  concerning  interference  problems  on  improved  networks 
need  to  be  validated  theoretically.  The  data  of  Mack  and  Robertson  or  of 
von  Seggern  cannot  be  used  for  prediction  because  they  represent  particular 
networks  of  stations  having  specific  noise  levels  and  also  because  they  are 
formed  In  response  to  necessarily  incomplete  epicenter  lists  created  with 
S ort  period  data  and  thus  cannot  be  Reneraiized.  Therefore  a simulation 
approach  is  required  to  fully  investigate  interference  effects  with  a variety 
o network  configurations  and  with  stations  having  lower  thresholds. 
order  to  perform  such  simulation,  knowledge  of  the  duration  and  shape  of 
typical  long-period  signal  trains  is  necessary.  Such  quantities  depend  on 
recor  i„g  distance,  magnitude,  depth  of  the  event,  and  structure  of  the  path 
we  can  parameterize  the  signal  according  to  these  variables,  then  a 
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reasonable  simulation  of  the  real  seismic  interference  situation  on  the  globe 
can  be  programmed.  Alsup  and  Lambert  (1973)  made  a statistical  characteriza- 
tion of  LR  (fundamental-mode  Rayleigh  waves)  signal  amplitudes;  but  because 
their  data  was  insuff  'ient  as  a global  sample,  because  they  only  measured 
signal  amplitudes  at  intervals  beyond  the  maximum  recorded  amplitude  (that 
part  of  the  signal  which  is  normally  termed  "coda"),  and  because  they  did 
not  take  a sufficiently  detailed  approach  to  signal  characterization,  we 
undertook  to  newly  measure  hundreds  of  long-period  signals  and  correct  these 
deficiencies.  The  first  part  of  this  report  describes  the  measuremnt  of 
these  signals,  the  second  part  describes  in  detail  the  character  of  LR 
signals  as  determined  by  these  measurements,  and  the  third  part  presents  the 
composition  and  results  of  a program  which  simulates  interference  on  a network 
of  long-period  stations.  We  will  be  concerned  with  Rayleigh  waves  (and 
therefore  vertical-component  recording)  only  in  this  report;  a parallel  study 
on  Love  waves  would  be  a more  difficult  task  since  horizontal  recordings 
always  contain  some  Rayleigh-wave  motion  which  would  need  to  be  eliminated. 
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SIGNAL  MEASUREMENTS 


High-gain,  long-period  recordings  of  the  LPE  (Long  Period  Experimental) 
network  were  used  for  signal  measurements  in  this  study.  The  particular 
sites  were:  ALQ,  CHG,  CTA,  EIL,  KIP,  KON,  OCD,  and  TLO.  Since  the  LPE 

instruments  have  high-gain  outputs  peaked  at  greater  periods  than  most  con- 
ventional long-period  instruments,  the  signal  shapes  reported  in  this  study 
should  differ  from  those  derived  from  conventional  instruments.  We  in  fact 
have  measured  several  signals  on  the  LPE  low-gain  seismograms,  which  were 
produced  by  instruments  with  responses  similar  to  those  of  the  WWSSN,  and 
found  that  the  signals  differed  distinctly  in  most  cases  from  those  measured 
on  the  high-gain  seismograms.  Thus  the  results  of  this  study  are  strictly 
applicable  to  LPE  high-gain  recording  only. 

In  order  to  sample  various  propagation  paths  on  the  globe,  we  selected 
eight  seismic  regions  (Flinn  and  Engdahl,  1965)  for  signal  sources;  Alaska- 
Aleutian  Arc  (1),  Mexico-Guatemala  Area  (5),  Andean  South  America  (8), 
Kermadec-Tonga-Samoa  Area  (12),  New  Guinea  (15),  Japan-Kuriles-Kamchatka  (19), 
Middle  East-Crimea- Balkans  (30),  and  Hindu  Kush-Pamir  (48).  Signals  actually 
measured  were  fairly  evenly  distributed  over  these  regions  and  over  the  eight 
stations  used  as  listed  in  Table  I.  A total  of  611  signals  made  up  the  final 
data  base.  Epicenters  in  1972  were  used,  and  only  those  events  whose  depths 
of  focus  were  listed  at  less  than  100  km  contributed  to  the  final  data  base. 
(Initially  some  deep  event  signals  were  measured,  but  their  shapes  were  seen 
to  differ  significantly  from  those  of  shallower  events.  Later,  a justifica- 
tion for  ignoring  deep  events  will  be  given,) 

We  were  interested  not  merely  in  that  part  of  the  signal  normally  termed 
"coda,"  that  is,  the  continued  motion  af^er  the  arrival  of  the  fundamental 
Rayleigh  mode  along  the  great  circle  path  but  also  in  the  relative  amplitudes 
of  the  entire  long-period  signal  train.  Our  objective  was  to  measure  the 
long-period  signal  from  the  P-wave  emergence  through  to  that  point  where  the 
coda  merged  into  the  background  noise.  For  this  reason,  data  selection  began 
with  the  largest  events  and  moved  down  into  smaller  events  only  if  they  were 
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TABLE  I 


r 


Distribution  of  Signal  Measurements 
by  Station  and  Seismic  Region 


Seismic  Region 


Station 

1 

5 

8 

12 

16 

19 

30 

48 

ALQ 

8 

10 

11 

9 

11 

10 

12 

9 

CHG 

12 

4 

5 

11 

10 

11 

5 

8 

CTA 

12 

7 

14 

9 

7 

11 

3 

3 

EIL 

13 

4 

11 

10 

8 

11 

12 

12 

KIP 

12 

9 

12 

10 

n 

12 

5 

2 

KON 

10 

9 

9 

10 

10 

12 

11 

12 

OGD 

11 

12 

11 

9 

12 

10 

9 

7 

TLO 

11 

12 

12 

9 

11 

10 

11 

4 

Seismic  Region  Name 

1 Alaska-Aleutian  Area 

5 Mexico-Guatemala  Area 

8 Andean  South  America 

12  Kermadec-Tonga-Samoa  Area 

16  New  Guinea 

19  Japan-Kur i les-Kamchatka 

30  Middle  East-Crimea-Balkans 

48  Hindu-Kush  and  Pamir 
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required  to  obtain  a sufficient  sample.  Many  of  the  largest  signals  overdrove 
the  recordings  and  were  thus  unusable.  In  order  to  get  a good  sample  for 
each  station  from  each  seismic  region,  it  was  necessary  to  use  some  signals 
for  which  body  waves  were  below  the  noise  level  and  finally  some  signals 
which  were  of  limited  duration  even  in  the  Rayleigh  and  coda  parts.  Arrival 
times  for  P,  S,  and  40-second  LR  were  predicted  by  routine  procedures  in 
order  to  pinpoint  t •»  arrival  of  these  phases.  Finally,  because  most  selected 
signals  had  no  visible  body  waves,  the  measurements  of  the  signal  train 
between  P and  LR  onset  were  not  subsequently  used  in  signal  characterization; 
but  some  data  is  available  for  a further  study  if  needed.  (Later  we  will 
justify  ignoring  this  part  of  the  long-period  signal  for  simulation  purposes.) 

The  analysis  procedure  was  to  measure  the  amplitude  and  period  of  the 
signal  at  two-minute  intervals  following  the  onset  of  P,  S,  or  LR,  whichever 
was  first  visible.  In  each  case  it  was  noted  with  which  phase  measurements 
began.  This  required  from  five  to  fifty  or  more  such  measurements  per  signal. 
The  amplitude  in  each  two-minute  interval  was  taken  to  be  the  maximum  within 
the  interval,  and  the  period  was  taken  to  be  the  dominant  or  average  period. 
The  measurements  were  coded,  punched,  and  edited  to  become  the  data  base  for 
long-period  signal  characterization. 


SIGNAL  CHARACTER 


DATA  BASE 

In  characterizing  the  long-period  signals,  we  disregarded  the  body-wave 
portion  and  concentrated  on  that  por'ion  following  onset  of  LR  because  body 
waves  are  detectable  only  for  event.,  above  roughly  4. 5-5.0 (n^)  and  because 
their  duration  is  short  relative  to  the  LR  and  coda  portions  on  the  long- 
period  seismogram.  The  portion  of  a typical  LPE  seismogram  taken  up  with 
body  waves  is  very  small,  only  a few  percent  of  a day's  recording  time.  Thus, 
any  precision  gained  by  attempting  to  model  the  body-wave  portion  of  the 
seismogram  for  a simulation  experiment  would  probably  be  small  and  would  not 
be  worth  the  added  complexity  introduced  into  the  simulation.  We  will  not 
be  significantly  underestimating  interference  effects  by  neglecting  body 
waves . 

We  first  plotted  the  signal  amplitudes  of  all  the  selected  events  for 
each  station-region  pair;  these  plots  are  contained  in  the  appendix.  Each 
event  was  normalized  by  the  maximum  amplitude  before  plotting;  note  that  the 
vertical  amplitude  scale  is  logarithmic.  These  plots  represent  the  signal 
from  LR  onset  to  whenever  the  coda  died  out,  that  is,  became  indistinguishable 
from  background  noise.  In  most  cases  the  plots  begin  with  the  longest  period 
LR,  which  is  normally  near  40  sec  or  greater  for  large  events  recorded  by 
these  LPE  instruments.  The  maximum  signal  amplitude  generally  occurs  within 
3 or  4 minutes  of  LR  onset  and  occurs  after  10  minutes  in  only  a few  cases. 
Some  signals  were  used  which  had  low  enough  S/N  ratios  that  the  earliest  LR 
cycles  of  motion  with  periods  of  40  seconds  or  more  were  probably  not  seen; 
in  these  cases  the  measurements  began  shortly  before  or  at  the  maximum 
amplitude.  Since  this  maximum  occurs  close  to  the  earliest  LR  arrival  anyway 
and  since  the  long  coda  train  will  not  be  shifted  too  much  relative  to  larger 
events,  it  is  acceptable  to  plot  these  smaller  events  with  the  first  measured 
arrival  at  zero  time  even  though  it  is  not  strictly  the  fundamental  LR  onset. 

If  all  events  for  a station-region  were  of  equally  high  S/N  ratio,  we  would 
expect  more  similarity  among  signal  shapes  for  each  station  region  pair.  We 
might  have  obtained  the  same  improvement  by  trying  to  align  signals  according 
to  reported  periods  before  plotting. 
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For  each  station-region  pair  in  the  appendix,  an  average  signal  shape 
is  computed  using  for  each  two-minute  interval  the  equally  weighted,  normalized 
amplitudes  trom  all  the  events  having  a measurement  at  that  interval.  The 
average  coda  shape  so  derived  is  less  reliable  as  time  increases  and  the 
measurements  become  fewer.  These  average  coda  shapes  are  reproduced  and 
grouped  by  station  in  Figure  1.  We  provide  in  Table  1.1  the  average  epicentral 
distance  for  events  in  each  station-region  pair  and  the  rough  percentage  of 
oceanic  and  continental  path  involved.  These  numbers  relate  to  the  coda 
shapes,  as  will  be  discussed  next. 

PARAMETER IZ 1 NG  SIGNAL  SHAPES 

We  next  attempted  to  parameterize  the  average  signal  shapes  shown  in 
Figure  1.  Only  nonlinear  functions  were  considered,  and  a program  which 
fitted  these  functions  in  a least  squares  sense  by  Gauss-Newton  iteration 
on  the  parameters  was  employed.  As  seen  in  the  plots,  every  coda  ends  with 
some  finite  amplitude,  presumably  the  background  noise  level;  therefore, 
before  fitting,  we  added  one  point,  equal  in  amplitude  to  this  last  point, 
to  the  beginning  of  each  set  of  amplitudes  to  be  fitted,  at  one  minute  prior 
to  the  first  measured  amplitude.  The  time  scale  as  seen  in  Figure  1 was  then 
shifted  back  one  minute  for  fitting  purposes.  This  results  in  a rather  rapid 
rise  time  to  the  signal.  Although  this  rapid  rise  will  not  realistically 
reflect  most  signal  shapes  because  body  waves  and  higher-mode  LR  usually 
precede  the  fundamental  LR,  we  feel  it  is  not  unrealistic  enough  to  degrade 
our  simulation  of  long— period  interference  since  the  amplitudes  preceding 
fundamental  LR  are  usually  much  smaller  and  in  many  cases  are  simply  background 
noise. 

Several  functional  forms  were  used  in  attempting  to  rit  the  normalized 
coda;  by  far  the  most  satisfactory  function  with  £ 4 parameters  was: 

P2C  P3C 

1°K10  IMt)]  = p^(l-e  ) e + P4 

where  t is  measured  in  minutes  and  log  A is,  arbitrarily,  the  logarithm  of 
the  average  of  the  normalized  codas  seen  in  the  appendix  plus  0.3.  This 
function  includes  an  exponential  decay  term  eP-^  with  scaling  parameter  p^; 
the  (l-eP2t)  term  provides  a smooth  rise  from  the  noise  background  amplitude 
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(RELATIVE  AMPLITUDE)  LOG..  (RELATIVE  AMPLITUDE) 


TIME,  min 


TIME,  min 


Figure  1.  Average  signal  shapes  (LR+coda)  from  shallow  epicenters 
for  64  station-region  pairs. 


-8- 


i 


T3 

Cj  u 
u — < 
3 C 


* ft* 
C3 

0)  c 


o 


u x 

o CJ 

^ y: 

I 

i*  c 

u c 
a — « 


U rj 

a *J 
u X 


X XT 


!3 

w ^ 


X 

<; 

h* 


C X 
^ C 


0)  X 
U C 


*-j  a 

« x 
o X 


<y  r: 

x c 
cj  x 

u 

x 

> 


o 

C 

O 

O' 

CN 

X 

\ 

\ 

\, 

o 

0 

o 

m 

v£ 

m 

\D 

v£> 

X 

-10- 


T 


— ** 


at  zero  time  to  the  maxinuim  amplitude;  and  the  additive  term  p^  relates 
directly  to  the  logarithm  of  the  background  l.oise  amplitude  and  does  not 
pertain  to  coda  shape  itsell.  The  final  values  o!  the  four  parameters  for 
each  itatiou-regiun  pair  are  given  in  Table  III  along  with  the  residual  mean 
square  error  between  the  estimated  l'unciton  and  the  actual  data.  (These 
particular  parameters  represent  only  the  normalized  signals,  which  will  be 
appropriately  scaled  when  simulating  actual  seismic  observations.)  Seven 
cases  where  convergence  did  not  occur  are  indicated  in  the  table;  however, 
the  residual  square  error  in  these  cases  is  apparently  not  any  larger  than 
lor  the  57  converging  cases,  and  so  the  tit  is  considered  good  for  even  these 
nonconverging  cases.  Several  cases  are  also  noted  where  p^  hit  the  boundary 
limit  and  was  held  constant  through  subsequent  iterations;  again,  the  resi- 
dual error  ir  those  cases  is  typically  no  greater  than  those  cases  where  p^ 
was  allowed  to  vary  throughout. 

It  was  thought  that  the  parameters  p^,  and  p^  should  show  some 

dependence  on  the  epicentral  distance  at  which  signals  were  recorded.  The 
average  distance  for  each  station-region  pair  was  given  in  Table  11.  In 
figure  2 p^,  p^,  and  p^  are  plotted  versus  these  distances.  A slight  depen- 
dence on  distance  is  apparent  for  each  parameter,  and  conjectural  lines  have 
been  visually  fit  to  the  data.  In  each  case  the  dependence  is  as  would  be 
expected  due  to  the  effects  ol  dispersion  on  Rayleigh  waves;  that  is,  the 
fundamental  I.K  is  progressively  stretched  out  in  time  with  increasing  epicentral 
distance,  leading  to  slower  rise  and  decay  times  for  the  envelope  of  the  signal. 
These  parameters  displayed  in  Figures  2a,  2b,  and  2c  are  calculated  for  the  coda 
together  with  the  fundamental  LK  arrival;  and  so  a large  scatter  with  distance, 
unrelated  to  dispersion  characteristics,  occurs.  (liven  the  data  in  these 
figures,  the  use  of  :i  constant  p^,  p , or  p^,  say  the  90°  value,  would  not 
lead  to  a gross  misrepresentation  of  signal  shape  at  any  distance  since 
scatter  at  the  extremes  of  0°  and  180°  seem  to  overlap  in  all  three  cases. 

DURATION 

For  simulation  purposes  it  is  necessary  to  estimate  an  appropriate 
duration  for  the  signal,  related  to  tiie  event  magnitude  and  the  epicentral 
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TABLE  111  (Continued) 
Signal  Parameter  Estimates 
P,1  P,t 
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SIGNAL  DURATION  (min)  SIGNAL  DURATION  (min) 


0 20  40  60  80  100  120  140  160  180 


Figure  Ab.  Duration  of  signals  where  1.00  _ log  A/T  < 1.30  (A/T  is 
in  units  of  millimicrons/seconds)  versus  the  epicentral  distance  of 
the  recording. 
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LONG-TERM  CODA  AMPLITUDES 


As  stated  previously,  many  signals  could  not  be  used  for  measuring  signal 
shape  because  they  were  overdriven  on  the  LPE  recordings.  With  this  restric- 
tion we  were  unable  to  measure  the  long-term  coda  amplitudes  visible  for 
large  events  our  longest  signal  in  the  data  base  was  132  minutes  while  codas 
for  large  events  (m^  >5.5  approximately)  continue  for  several  hours.  Although 
events  with  > 5.5  occur  at  the  rate  of  only  roughly  one  every  four  days, 
their  prolonged  coda  makes  them  account  for  a much  more  significant  portion 
of  the  interference  problem  than  their  mere  number  might  indicate.  We  felt 
it  necessary  to  gather  some  data  on  the  coda  amplitudes  for  large  events  to 
be  used  whenever  needed  in  the  simulation  program;  a simple  extrapolation  of 
the  signal  shape  function  estimated  from  smaller,  shorter  signals  might  be 
very  inaccurate  for  times  greater  than  an  hour  or  so,  especially  since  we 
have  used  an  exponential  function  in  the  fitting  of  signal  shapes. 

To  this  end  we  examined  our  LPF,  eight-month  data  file  (von  Seggern,  1974) 
for  LR  signals  which  had  been  indicated  as  overdriven  at  a station  by  the 
analysts.  Signals  were  then  examined  to  see  both  if  they  extended  for  at 
least  200  minutes  after  the  LR  arrival  and  if  they  were  readable  within  10 
to  30  minutes  after  this  arrival.  Also,  it  was  necessary  that  no  signals 
from  other  events  large  enough  to  be  recorded  at  the  station  arrived  during 
this  coda.  Clearly,  subjective  judgement  was  involved  here,  but  we  tried 
to  be  reasonably  sure  that  no  other  events  were  contributing  to  the  coda. 

These  criteria  together  were  rather  stringent;  seventeen  signals  which  met 
them  were  found  though,  and  they  were  from  six  events  as  listed  in  Table  IV. 
Amplitude  and  period  of  the  LR  signal  and  coda  was  measured  at  ten-minute 
intervals  after  the  LR  onset,  starting  whenever  the  recording  became  readable, 
in  the  same  manner  as  done  at  two-minute  intervals  for  the  611  signals  in  our 
main  data  base  (except  that  amplitudes  associated  with  periods  > 50  sec  were 
ignored  here). 

Recorded  amplitude  (millimeters)  versus  time  for  the  six  events  is  shown 
in  Figure  5.  Arrivals  of  R2  are  apparent  in  many  of  the  plots  where  amplitudes 
increase  with  time,  and  these  arrivals  could  indeed  be  seen  on  the  recordings 
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Epicenter  Information  (NOS)  on  Events 
Used  for  Long-Term  Coda  Measurements 
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Figure  5.  Amplitude  versus  time  after  LR  onset  for  large  earthquakes 
(LPE  vertical  component). 


Figure  5.  Amplitude  ver 
(LPE  vertical  component) 


wheie  usually  longer  periods  accompanied  the  amplitudes.  Using  only 
amplitudes  at  _30  minutes  and  grouping  all  signals  for  all  events  together, 
a computed  regression  of  log  ^ 011  in  minutes  had  a slope  of  -.006  7. 

This  line  has  been  drawn  on  all  six  figures,  and  it  clearly  fits  all  the 
coda  decays  at  all  times  to  within  a factor  of  approximately  two  if  it  were 
shifted  up  or  down  by  an  amount  appropriate  to  each  of  the  six  signals. 

Using  the  decay  computed  for  large  events'  Coda  here  and  the  signal  shape 
calculated  from  smaller  events'  signals,  we  are  able  now  to  generate  signal 
shapes  from  LK  onset  out  to  several  hours  if  necessary  for  our  simulation  of 
interference  problems. 

CODA  PERIODS 

The  frequency  of  coda  waves  is  important  since  detection  of  newly- 
arriving  signals  depends  not  only  on  their  amplitude  relative  to  the  coda 
buc  also  their  frequency  content  relative  to  the  coda.  Mosf;  of  the  LR  signals 
detected  on  the  LPE  network  had  a visible  bandwidth  at  least  from  .025  to 
.05  Hz  (von  Seggern,  1974).  Due  to  dispersion,  the  energy  in  this  bandwidth 
has  a characteristic  time  distribution  which  further  enhances  the  detect- 
ability of  LR.  Since  our  signal  analysis  for  this  study  included  reporting 
of  the  period  at  each  two-minute  interval,  we  are  able  to  present  the 
relation  of  dominant  period  to  the  coda  time.  Figure  6 shows  plots  for 
typical  recordings  with  high  S/N  ratio.  Again,  the  zero  time  is  the  funda- 
mental LF  onset.  Note  that  each  of  the  eight  stations  is  represented  and 
that  each  of  the  eight  seismic  regions  is  also;  so  there  should  be  no 
particular  bias  to  this  sample.  The  salient  feature  of  these  plots  is  that 
except  for  a few  minutes  during  the  arrival  of  the  fundamental-mode  LR  at 
the  start  and  perhaps  the  arrival  of  this  phase  again  from  the  longer  great 
circle  path,  the  dominant  period  of  the  signals  is  quite  limited,  centered 
at  roughly  20  seconds.  This  coda  characteristic  implies  that  newly-arriving 
signals  from  other  events  should  be  easily  detectable  visually  if  their 
amplitude  is  at  least  as  large  as  that  of  the  coda  because  of  their  band- 
width and  dispersion.  Moreover,  if  processing  is  feasible,  a narrow  norch 
filter  or  a low-pass  filter  with  corner  at  .05  Hz  may  be  helpful  in  extracting 
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Figure  f>.  Typical  pints  of  dominant  signal  period  versus  time  alter 
onset  of  fundamental-mode  I.K.  •,< 
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Figure  6.  Typical  plots  of  dominant  signal  period  versus  time  after 
onset  of  fundamental-mode  LR. 
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signals  from  rhe  coda  of  other  events.  Because  at  least  75%  of  the  time 
span  of  the  signals  is  dominated  by  this  .05  Hz  coda  as  shown  in  Figure  5, 
there  is  cause  for  optimism  in  recovering  by  digital  or  analog  processing 
many  events  which  are  visually  masked  within  the  coda  of  larger  events  at 
single-sensor  sites. 

Further  evidence  for  this  20-second  dominance  in  coda  is  obtained  from 
our  17  signals  with  > 200-minute  durations.  Out  of  all  the  measured  periods 
on  these  signals,  three- fourths  were  in  the  range  of  18  to  25  seconds.  There 
was  no  discernable  systematic  change  in  period  over  the  duration  of  these 
long  codas,  aside  from  the  lengthening  of  period  often  accompanying  arrival 
of  f^.  (We  note  that  codas  on  the  low— gain,  WWSSN— type  responses  had 
predominate  periods  somewhat  shorter,  15  to  20  seconds,  than  the  codas  on 
the  high-gain  counterparts.) 

The  fact  that  the  dominant  periods  in  the  coda  are  mostly  near  20  seconds 
allows  us  to  infer  an  interesting  geophysical  datum,  Q,  for  this  period. 
Remembering  that  we  found  log^  A u — .0067t  for  coda  decay,  taking  the 
antilog  and  converting  to  base  e,  and  then  assuming  an  average  3.2  sec/km 
group  velocity  for  20-second  LR  waves  over  the  glove,  we  obtain 
A " exp(-.000080r)  where  r is  ;n  kilometers.  Inelastic  attenuation  of  surface 
waves  is  commonly  expressed  as  A « exp(-^r/QUT)  where  Q is  the  quality  factor, 
L is  group  velocity,  and  T is  period.  Now,  equating  exponents  and  substi- 
tuting in  our  values  for  U and  T,  we  get  Q : 630.  Note  that  we  have  not 

included  the  effect  of  dispersion  on  amplitudes,  a factor  proportional  to 

,-1/2  J 

, and  so  the  present  estimate  of  Q is  conservatively  low.  Other  effects 
must  be  taken  into  account;  for,  clearly,  20-second  waves  late  in  the  coda 
must  have  travelled  extremely  complex,  circuitous  paths  (albeit,  excellent 
ones  for  LR  propagation)  wherein,  to  be  precise,  geometrical  spreeding,  mode 
conversion,  and  reflection-refraction  coefficients  would  need  to  be  calculated 
for  any  theoretical  model.  But  since  the  arrivals  in  the  coda  have  these 
diverse,  unknown  paths,  such  effects  would  be  very  difficult  to  quantify; 
we  know  only  the  arrival  time  of  each  amplitude  measurement  and  thus  the 
length  of  the  path  travelled. 
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SIMULA  I ION  MODEL 


We  desire  to  simulate  in  the  time  domain  the  seismic  activity  seen  by 
typical  LPE  sites  and  thus  by  the  network  as  a whole.  Some  events  will 
produce  signals  above  the  noise  level  while  others  will  not,  according  to 
signal  and  noise  amplitude  levels  at  each  particular  site.  For  those  signals 
which  are  above  the  noise,  clear  detection  is  not  guaranteed  because  of  the 
possibility  of  signals  from  other  events  arriving  simultaneously.  The 
purpose  of  our  simulation  is  to  determine  what  percentage  ot  events  at  a 
given  magnitude  are  thus  affected  when  we  vary  network  threshold  and  con- 
figuration. Empirical  measures  have  been  reported  by  von  Seggern  (1974)  for 
the  LPE  network,  and  these  empirical  results  will  aid  in  calibrating  our 
simulation  program;  that  is,  we  will  attempt  to  reproduce  those  results, 
which  grew  from  visual  analysis  of  seismograms,  with  a simulation  model 'run 
on  the  computer  before  proceeding  with  other  cases. 

In  the  calibration  stage  it  is  necessary  to  run  the  model  with  a real 
epicenter  list,  in  fact  the  very  list  with  which  the  LPE  data  base  was 
formed.  But  if  we  are  to  employ  the  model  for  very  sensitive  networks  with 
thresholds  lower  than  achieved  so  far  in  seismology,  we  must  simulate  all 
events  that  could  be  detected  with  even  a very  small  probability  since  each 
event,  no  matter  how  small,  can  possibly  contribute  to  the  interference 
problem  if  it  is  near  a station.  No  presently  available  epicenter  list  has 
a low  enough  threshold  to  assure  that  all  events  which  might  produce 
observable  long-period  recordings  are  entered.  (Recordings  from  the  LPE 
sites  typically  show  five  or  more  events  per  day  than  are  on  the  NOS  list.) 
Therefore,  in  the  full  simulation,  seismic  events  will  be  generated  by  a 
Monte  Carlo  technique.  In  this  way,  although  we  are  able  to  create  the 
necessary  activity  unreported  on  real  epicenter  lists,  we  replace  the  known 
pattern  of  earthquake  occurrence  in  the  earth  with  a random  (except  for  the 
existence  of  aftershock  sequences)  spatial  and  temporal  variation;  but  this 

departure  from  reality  should  not  significantly  detract  from  the  utility  of 
the  results. 
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magnitude  > Mg.  Mg  is  the  proper  magnitude  to  simulate  by  Monte  Carlo  since 

it  relates  to  signal  amplitudes  as  measured  for  this  study.  No  events  with 

Ms  < 3,0  are  considered  since  only  rarely  are  surface  waves  from  an  event 

below  3.0  recorded  beyond  local  or  regional  distances,  where  only  a few 

cycles  of  signal  are  visible  above  the  background  anyway.  Clearly,  events 

of  Ms  < 3.0  cannot  cause  significant  interference.  The  above  relation  implies 

an  average  waiting  time  of  27  minutes  between  events  with  M > 3#0. 

s — 

It  is  very  important  to  include  aftershocks  in  our  earthquake  simulation 
since  aftershocks  will  produce  more  mixing-masking  incidents  than  random 
main  events.  Shlien  and  Toksoz  (1970)  have  shown  that  earthquake  lists  can 
be  described  by  a generalized  Poisson  process  wherein  a random  number  of 
aftershocks  vmost  often  zero)  is  associated  with  each  main  shock.  This  ran- 
dom number  of  aftershocks  is  selected  from  the  discrete  Z (or  Pareto)  density: 

P (n)  = c”1(E)t."E 

where  ;(E)  is  the  Riemann  zeta  function.  For  n = 1,  there  is  just  a main 
event  and  no  aftershocks.  Calculations  based  on  the  above  density  with 
E = 3.0,  a value  found  for  NOS  epicenters  from  1963  to  1968  by  Shlien  and 
Tokstiz,  predicts  that  .23  of  all  earthquakes  will  be  of  aftershock  nature  and 
.77  will  be  main  shocks.  Sadeh  and  Meidav  (1973),  applying  very  liberal 
aftershock  criteria  to  NOS  epicenters  from  1964  through  1972,  found  that 
roughly  17%  of  the  events  were  aftershocks  so  that  the  true  percentage  of 
aftershocks  is  probably  somewhat  less.  Since  the  seismicity-magnitude 
relation  given  above  implies  an  average  waiting  time  (WT)  of  27  minutes 
between  origin  times  for  all  events  with  Mg  > 3.0,  the  waiting  time  for 
main  shocks  only  must  be  revised  upward.  The  relation 

27  < .77  (WTMS)  + .23  (WT^) 

must  hold  for  the  average  waiting  times  of  main  shocks  (WT  ) and  after- 
shocks (WTas).  The  equality  would  hold  only  if  WT^  approached  zero. 

Aftershock  waiting  times  are  known  to  vary  both  with  magnitude  of  the  main 
shock  and  with  time  after  the  main  shock  (Shlien  and  Toksoz,  1974);  but 
there  is  no  analysis  of  aftershock  waiting  times  on  general  basis  known  to 
the  author.  Consequently,  we  arbitrarily  assumed  WT^  to  be  15  minutes, 
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tliu.s  n‘lt"11'  “r|>&  - 11  (In  actual  simulations,  wr  - 32  gUVe  t|lt, 

.orro.t  average  vail  I an  t imo  „f  27  minutes  for  nil  events.')'  Because  we 

ronard  this  average  waiting  time  lor  aftershocks  to  he  on  the  low  side  and 

because  the  number  of  alters!,, .cits  seems  to  be  too  hiBh  in  the  model  of 

Toksox  and  Shlicn  <1970).  we  expect  to  certainly  not  underestimate  the  amoun, 

o,  inter, erence  on  a network  when  wc  incorporate  aftershocks  into  our  simu- 
la t ion. 

Assumint  that  both  main  events  and  aftershock  sequences  associated  with 
a particular  main  event  occur  as  Poisson  processes,  their  waiting  times  are 
exponentially  distributed  and  can  bo  randomly  generated  by, 

m " CTMS  ‘ n(1/*>  or  KT  - WTas  . !:n(l/x) 

where  x is  a random  variable  drawn  from  a uniform  distribution  (o-x-1).  (If 
no  aftershock  sequences  are  generated,  we  revert  to  an  average  waiting  time 
Of  27  minutes  as  stated  above.)  Assuming  to  be  distributed  according  to 
the  above  sesimlcity-magnitude  relation,  a random  M (>3.0)  can  be  generated 
for  main  events  and  aftershocks  by: 

Ms  = I n(l/x)/ (.872  • ;.n  10)]  + 3.0 

whore  x again  is  randomly  drawn  from  the  uniform  distribution  (0x1)  For 
altershocks,  we  do  not  allow  a randomly  generated  to  exceed  that  of  the 

**“  ■',n‘l  " “*•  draw  new  random  x’s  until  the  required 

number  of  altershocks,  n-1,  is  attained  with  this  constraint.  This  procedure 
Wtll  make  aftershock  sequences  for  small  events  appear  to  resemble  swarms 
since  we  already  require  that  all  events  have  * > 3.0.  For  location  of  the 

mxnn  event,  random  longitude  and  latitude  are  generated  by: 

J.ON  = 360x  - 180 
_ 180 

- — — arcs  in  x 

Where  each  x is  a different  random  number  drown  from  the  uniform  distribution 

Another  random  number  x is  drawn,  and  if  x > .5,  the  sign  of  LAT  is  changed, 

thus  producing  an  equal  number  of  events  in  the  northern  and  southern 
hem i spheres . 
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When  a real  epicenter  list  is  used  in  providing  sources  for  the  simula- 
tion run,  a conversion  from  mfa  to  ^ must  he  made  for  each  event  since  events 
are  touti)iely  assigned  but  only  occasionally  Mr.  We  used 

M„  = 1.45  m.  - 2. 34 
*>  b 

a regression  counted  by  von  Soggern  (1974)  using  over  3000  M -m  poilus  fm. 
global  earthquakes,  s 3 

We  have  previously  discussed  our  reasons  for  neglecting  the  body-wave 
portion  of  the  seismogram  in  simulating  long-period  signals-relative  to 
surface  waves,  they  simply  do  not  occupy  much  time  on  typical  recordings. 

**  with  neglecting  deep  events  and  small  events,  this  will  tend  to  cause  an 
underestimation  of  interference  effects.  As  mentioned  previously,  these 

omissions  should  not  seriously  bias  the  results  of  the  simulation  program 
though. 

DETECTION  WITHOUT  INTERFERENCE 

By  generating  spatially  and  temporally  random  events  and  assigning  them 
signal  amplitudes  and  shapes  at  calculated  epicentral  distances,  we  are  able 
to  survey  a station  or  network  in  the  simulation  program  for  its  response  in 
terms  of  detections,  non-detections,  and  interference  cases.  In  the  case  of 
-signal  interference,  the  shapes  will  be  important;  but  let  us  first  consider 

the  simple  response  where  a signal  arrives  alone  and  only  amplitude  need  be 
considered  for  detection. 

Each  station  used  in  the  simulation  run  is  assigned  a background  noise 
value.  This  value  is  the  median  of  the  maximum  zero-to-peak  noise  amplitudes 
(millimicrons)  measured  in  numerous  40-minute  windows  of  seismic  background 
noise  (von  Seggern,  1974)  if  an  actual  I.PE  site  is  used  or  estimated  value 
(50m.)  if  a hypothetical  site  is  used.  No  noise  variance  is  allowed;  that 
is,  the  background  is  assumed  to  remain  constant.  Signal  amplitudes  at  each 
site  are  calculated  by: 

M -1.06  logA  - 1.43 

A = T • 10 

Where  A ls  a *er.°-to-penk  amplitude  (millimicrons),  T is  assumed  to  be  a 
constant  20  seconds  period,  and  £,  is  distance  in  degrees.  (The  M formula 

R 
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implied  here  is  the  result  of  a study  of  LPE  data  in  von  Seggern,  1974.)  No 
signal  variance  is  allowed  either;  so  we  do  not  attempt  to  model  recognised 
magnitude  variation  due  to  radiation  pattern  and  path  effects.  A signal  is 
declared  as  detected  at  a station  if  its  amplitude  as  calculated  above 
exceeds  the  background  noise  amplitude.  Thus  detection  jumps  from  zero  to 
100  percent  at  S/N  ratio  ol  1.0,  a value  determined  by  von  Seggern  (1974)  for 
Rayleigh  waves  on  U’E  recordings. 

Our  restriction  of  zero  noise  and  signal  variance  requires  some  dis- 
cussion. In  a reel  detection  situation,  both  signal  and  noise  variance 
contribute  to  the  cumulative  normal  shape  of  observed  detection  curves. 

Since  the  rate  of  seismicity  increases  with  decreasing  magnitude,  such  real 
detection  curves  entail  many  detections  at  magnitudes  lower  than  the  cutoff 
on  the  artificial  unit-step  detection  curve  which  we  have  assumed  for  our 
simulation.  Clearly  our  unrealistic  assumption  would  adversely  affect  the 
simulation  results  for  detection  and  interference  if  seismic  activity  were 
confined  to  one  point  on  the  globe.  However,  we  in  fact  simulate  random 
activity  over  the  entire  globe,  and  large  signal  deviations  for  a given  Mg 
occur  naturally  at  all  stations  because  of  the  large  scatter  in  distances 
to  events  with  that  M . This  scatter  tends  to  overwhelm  that  which  would 
result  if  we  incorporated  signal  variance  for  events  and  noise  variance  at 
stations.  We  therefore  conclude  that  our  restriction  will  not  seriously 
bias  the  results. 

DETECTION  WITH  INTERFERENCE 

If  two  or  mare  signals  interfere,  we  still  check  each  signal  for  detec- 
tion as  in  the  simple  non-interference  case;  however,  if  each  is  above  the 
noise  level,  some  additional  logic  must  be  employed  to  decide  whether  they 
should  be  classified  as  "mixed"  or  "masked"  rather  than  "detected."  In  a 
real  situation  many  events  which  would  be  detected  when  arriving  alone  must 
be  relegated  to  the  mixed  or  masked  category  by  the  analyst.  The  simulation 
program  must  emulate  the  analyst's  logic  as  closely  as  possible.  This  is  a 
formidable  problem  since  the  situations  in  which  events  interfere  form  an 
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extensive  and  varied  fabric  to  which  the  analyst  applies  visual  criteria, 
often  vague  and  intuitive,  in  deciding  whether  a particular  event  is  seen  or 
not.  This  contrasts  with  the  fairly  straightforward  decision  process  when 
an  event  arrives  singly.  The  simulation  program  must  apply  algorithms  with 
quantified  decision  levels  to  signals  which  are  also  well  quantified.  Clearly, 
it  is  in  this  case  now  that  signal  shape  will  be  important. 

The  dependence  of  signal  shape  on  epicentral  distance  was  previously 
presented.  We  concluded  that  this  dependence  was  rather  weak;  nevertheless 
in  the  simulation,  distance  dependence  is  crudely  modelled  in  that  a signal 
is  assigned  parameters  p^,  p^  and  p^  equal  to  values  at  15°,  60°,  and  135° 
on  the  straight  lines  in  Figures  2a,  2b,  and  2c,  depending  on  which  distance 
range  (0-30  , 30-90°,  or  90-180°)  the  station  fell  in.  Note  that  for  the 
Pi  values  there  is  a large  scatter  which  we  do  not  choose  to  model;  we  will 
return  to  this  matter  later.  Since  we  often  need  to  generate  signals  with 
durations  longer  than  those  measured  for  signal  shape  and  since  we  have 
shown  that  the  latter  part  of  signal  codas  can  be  reasonably  approximated  by 
an  exponential  decay,  we  model  that  part  of  the  coda  beyond  30  minutes  after 
LR  onset  by  the  decay  relation  computed  for  the  17  large  magnitude  events: 

log10A  " log10Ao  = ■•°067  (t“ 30) 

where  Aq  is  the  amplitude  calculated  at  30  minutes  and  t is  the  time  in 
minutes  after  LR  onset.  Since  the  maximum  of  the  signal  envelope  is  tied 
to  Mg  and  thus  the  calculated  maximum  amplitude  at  the  station,  amplitude 
in  millimicrons  at  any  point  in  the  signal  can  easily  be  calculated  by  scaling 
the  peak  of  the  signal  envelope  to  this  maximum  amplitude. 

A check  on  the  validity  of  our  signal  shapes  in  the  simulation  model 
can  be  made  by  comparing  the  durations  of  simulated  signals  with  real  ones. 
Figure  3 presented  signal  duration  versus  maximum  amplitude  for  the  real 
events  that  were  measured  for  signal  shape.  A preliminary  run  of  the  simula- 
tion program  was  made,  with  signals  modelled  as  described  above,  to  produce 
a similar  body  of  data,  shown  in  Figure  7.  The  straight-line  visual  fit  to 
the  data  of  Figure  3 is  reproduced  here,  and  also  the  general  extent  of  the 
Figure  3 data  is  shewn  by  a dashed  perimeter.  The  simulated  durations  (that 
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Equating  to  zero  and  solving  for  t gives 

m ° 


Using  the  values  of  d at  1 1°  f.n° 
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thC  V:*lues  >1  «"•  lm  <*>«*..  results  In  the  portions  „f  t,„ 

curves  t„  Figure  7 up  to  thirty  .!«««.•  duration.  For  lo.,Ber  hurutions. 
the  long-term  coda  relation  estahlished  earlier  must  he  used  In  con), .notion 
With  Che  short-term  relation.  11, o long-term  relation  is  again: 
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Multiplying  this  hy  the  l-term  ratio  given  above  and  neglecting 


results  in: 
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l.smg  the  same  values  ot  p.  and  Lm  as  above  produces  the  portions  of  the 
three  curves  in  Figure  7 for  durations  greater  than  thirty  seconds.  The 
actual  durations  "observed"  from  the  simulation  run  lie  outside  these  curves 
in  some  instances  simply  because  the  noise  level  at  each  station  was  different 
' r°m  50  zero-to-peak.  Ove.all,  the  calculated  duration- amplitude  rela- 
tions in  Figure  7 do  not  violate  the  real  data  of  Figure  3,  and  no  real  data 
was  available  to  coni  inn  our  prediction  of  a reduction  in  the  slope  of  the 
duration-amplitude  relation  for  very  large  events  with  maximum  A/T  > 1000. 

The  durations  in  Figure  7 do  not,  as  mentioned  previously,  reflect  the 
obseived  scatter  in  durations  of  actual  signals.  This  is  due  to  the  con- 
straints of  the  p.  which  we  have  made  for  the  purpose  of  limiting  computational 
time  ol  the  simulation  program.  With  the  three  sets  of  p.'s  used  for  three 
distance  ranges,  the  program  can  preset  the  three  corresponding  signal  shapes 

and  merely  scale  them  up  or  down  according  to  the  current  event  M ; but  if 

s 
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Zulu  7 KU~*d  “ 7"*  “atl*tl“1  «•»  signal  shape  for 

simnl  ti"J  J'"V'  U"  "1,iat  1“V“1  l°  bL'  c™l',,,-cd-  “«  did  make  one 

elation  run,  hoover,  fo  test  the  effect  of  varying  the  p.’s  to  more 

cosely  model  actual  variability  in  signal  shape.  this  run  each  signal 

7 /0i“  U'VCl  “ “*«-  - -!•  -V  - P3  imposed  of  a term  equal 

he  value  on  the  straight  lines  in  F1g„res  2a.  2b.  and  2c  at  the  proper 

-oce  plus  a random  term  reflecting  the  statistical  distribution  of 
po  nts  vertically  about  a point  on  the  line  in  these  figures.  This  distribu- 
tion was  assumed  to  be  normal  in  each  case,  with  the  mean  on  the  line  and 
with  standard  deviations  of  .45.  1.0,  and  .020  for  p p and  „ 

tively.  (No  positive  values  of  p p or  | * 2 * 3*  respec~ 

Pl’  P2*  or  P3  were  allowed  though;  they  were 

ma  negative  values  if  this  occurred.)  This  refined  model  was 

used  for  the  first  thirty  minutes  of  signal  only;  the  constant  -.0067 
exponential  decay  for  the  long-term  codas  was  retained  without  variation. 

.e  s mulation  tun  with  variable  signal  shapes  resulted  in  roughly  202  more 
xing  and  masking,  and  consequently  fewer  detections,  than  a run  identical 
y respect  except  that  signal  shapes  were  constrained  to  the  three 
sets  of  constant  p^s  for  three  distance  ranges.  However,  the  required 
computer  run  time  was  quadrupled.  We  do  not  feel  that  the  refined  modelling 
jus  ies  this  increase  in  expenditure,  and  the  results  given  in  this  report 
will  be  for  the  constrained  model  of  signal  shapes. 

«e  now  refer  to  the  logic  for  mixed  and  masked  events  presented  on  the 
simulation  program  flowchart  in  Figure  8b.  The  first  steps  of  the  program 
are  charted  in  Figure  8a  and  have  already  been  discussed.  At  this  point  we 
qualitatively  define  a "m.xed"  event  as  one  whose  signal  overlaps  the  signal 
Of  one  or  more  other  events  (considering  only  those  parts  of  the  signal  which 
exceed  the  background  noise)  and  which  is  detected  in  some  imperfect  sense 
such  that  measurement  of  period  and  amplitude  would  be  uncertain  due  to  the 
interference.  We  qualitatively  define  a "masked"  event  as  one  whose  signal 
exceeds  the  background  noise  for  some  time  Interval  but  which  is  effectively 
concealed  by  another  signal  arriving  at  that  time;  this  does  not  require  that 
the  amplitude  of  the  masking  signal  exceed  that  of  the  masked  signal  at  every 
point.  The  "mixed-masked"  logic  in  Figure  8b  is  an  attempt  to  emulate 
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Flowchart  of  interference  logic  in  simulation  program 


classification  of  mixed  and  masked  events  by  , h„m 

<»  logical  icc,  wore  'J*  " 

"“del  against  real  data  as  will  1 •*  i bating  tlio  sinulat  Ion 

*>c*.*. » logic  „ more  ^J|y  r ,he  ~*1~  "" 

occurrences  ore  covered  by  tho  , . " *"  ’•  s°tc  that  .ill  possibl, 

« - rectangular  dell , 17,  ^ - iagit 

^ ■ t22  • ’23-  *u  - *n>  =t  rr T"  *>  * ^ 

employ  curves  to  separate  Ll  23  ™0rC  3011  'Stlc  approach  uoohl 

separate  the  decision  areas,  but  we  resir  I ,h. 

os  not  only  indeterminable  but  also  unnecessary.  ‘"C"",U 


CALIBRATING  THE  SIMULATION  MODEL 


cop  .1:11::;  - i,ave  “ • — 

c^elr  coda.  We  have  checked  the  durl  i„7  ^Xus  R',ylei?h  “'1VeS  "" 
against  those  of  real  signals  w » generated  in  the  model 

•vv  - ompiitude-disri-re::::: 

i: ::  ~rp::;:,r::;:;7--s  - — 

-ok  to  x the  i:;::;  :;i;;87::r:riod  rec°rdi- »» « 

-on  program  to  that  reported  hy  analysts  working  wi7“  ^ ^ ^ 
Uml  «•»  *»  “"Posing  the  working  program.  J,  • “ ° 

factors  In  the  logic  scheme  of  Figure  9.  1"''°  Ve"  deteri"lnln«  th0 

for  our  calibration  data,  the  month  of  January  1972  „-s  , . 

flic  U'F.  data  lase  <v„„  Seggern.  1974,  ' «» 

detected,  mixed,  and  masked  for  stations'  77"  0"  m""bCr 
simulation  program  was  run  on  the  same  evil 

f ic ial , statistics.  The  parameters  T . g i " b“‘  “f“- 

varied  until  satisfactory  agreement  to  wit, 

, ’ L wit,lln  upproximately  20%  Wls 

attained  between  tb.  renl  and  simulated  results'  tbe  ■ , . 

»fked  and  masked  percentages  were  purpose!  t’  „ T 

Poccentages  since  deep  events  and  body  will  . ”'*7' 

simi.lation  The  mr  ' considered  in  the 

un’  lhe  parameter  values  ser  ,u ■ 

venues  set  in  this  manner  were: 
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= 6 minutes, 


T 

11 

= T 

12 

= T 

13 

T21 

= T22 

= T 

23 

A11 

= A13 

= . 25, 

A12 

= -5, 

A21 

= 2, 

A22 

= A23 

= 1. 

8 minutes, 


Admittedly,  this  is  not  a unique  set  of  parameter  values;  and  other  combina- 
tions of  significantly  different  values,  or  even  a change  in  the  basic 
decision  logic,  could  perhaps  provide  equally  satisfying  agreement  between 
the  simulated  and  empirical  data.  For  example,  decreasing  the  time  constant 
by  two  or  three  minutes  or  changing  the  amplitude  constants  by  502  does  not 
greatly  affect  the  results  of  the  simulation,  but  time  constants  of  1 minute 
or  amplitude  constants  twice  as  large  gave  unsatisfactory  results.  However, 

the  above  set  of  values  provides  satisfactory  agreement  and  has  the  property 
of  being  physically  reasonable. 


Before  proceeding,  we  also  need  to  cneck  the  seismicity-magnitude  graph 
of  the  simulated  data  against  the  relation  used  in  generating  M values  since 
the  addition  of  aftershock  sequences  may  significantly  distort  this  relation. 
In  Figure  10  we  present  sets  of  simulated  data  in  the  form  of  incremental 
graphs  of  number  of  events  versus  It  can  b,  seen  that  for  the  simulation 

with  no  aftershocks,  the  agreement  with  the  seismicity-magnitude  relation 
is  very  good  over  the  entire  range  while  for  the  simulation  including  after- 
shocks the  resulting  numbers  are  somewhat  too  high  at  the  lowest  magnitudes 
and  somewhat  too  low  at  the  intermediate  magnitude  due  to  a concentration  of 
aftershocks  toward  the  lower  magnitudes.  Even  for  the  simulation  with 
aftershocks,  the  variance  of  the  points  from  the  underlying  relation  is  not 
considered  to  be  significant  though,  and  a line  fit  through  these  points 
would  not  differ  much  from  the  given  one. 
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Seismicity-magnitude  plots  of 


simulated  earthquake  bulletins 


Figure  10. 


SIMULATION  RESULTS 


Using  values  as  previously  defined  for  all  the  parameters  in  the  program 
for  simulating  network  interference,  we  proceeded  to  obtain  results  for  hypo- 
thetical networks.  Five  such  arbitrary  networks  were  considered;  r>,  10,  and 
20-station  global  networks  and  5 and  10-station  northern-hemisphere  networks. 
The  -0  stations  ol  which  all  these  networks  are  comprised  ire  listed  in 
Table  V.  Note  that  each  network  covers  the  globe  (or  northern  hemisphere) 
somewhat  evenly.  Ten  of  the  locations  in  Table  V are  presently  LPE  sites; 
eight  of  the  remaining  ten  were  or  are  WWSSN  sites. 

fhe  simulation  results  for  the  five  networks  are  summarized  in  Table  VI. 
Here  station  results  are  averaged  over  all  the  stations  in  a particular 
network.  The  categories,  "detected",  "mixed",  and  "masked",  have  been 
defined  lor  a station,  but  not  as  yet  for  a network.  By  "detected"  on  a 
network,  we  here  mean  that  >_  two  stations  of  the  network  declared  a clear 
detection;  by  mixed"  that  ^ two  stations  declared  a detection,  clear  at 
none  or  at  most  one  station  and  mixed  at  the  remainder;  by  "masked"  that 
either  at  least  one  station  reported  the  event  as  masked  if  only  one  station 
reported  a detection  (mixed  or  clear)  or  at  least  two  stations  reported  the 
event  as  masked  if  no  station  reported  detection  (mixed  or  clear).  In  other 
words,  the  network  reports  an  event  as  mixed  if  signal  interference  hindered 
it  from  having  at  least  two  clear  station  detections,  and  it  reports  an 
event  as  masked  if  more  severe  interference  hindered  it  from  having  at  least 
two  detections  of  any  sort.  (We  mention  that  in  the  real  detection  situation, 
there  is  no  way  of  knowing  absolutely  whether  an  arrival  one  might  choose  to 
report  as  masked  would  have  been  detected  under  normal  background  conditions. 
The  simulation  program  is  artificial  in  this  sense  then:  it  knows  a pi iori 

v»ii.it  events,  when  interferred  with,  would  have  been  detected  without  inter- 
ference. On  the  other  hand,  the  real  data  is  imperfect  since  tile  "masked" 
label  will  be  arbitrarily  applied.  (The  number  of  reported  masked  events  in 
the  LPE  data  exceeded  the  number  of  truly  masked  events  because  the  analysts 
routinely  called  an  arrival  masked  if  its  m.  >_  4.0  and  A < 100°  and  if  a 
signal  from  another  event  was  present.  A crude  adjustment  to  the  empirical 
data  was  made  in  calibrating  the  simulation  program  by  noting  the  observed 
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tabu;  v 


>t.it  ion  Local  ions  t'sed  in  the  Simulation 

Station  Code 

Locat ion 

Lnt i tu  ie 

Lonir  i t tido 

AF  1 

At  iamalu,  S.im.io  Islands 

1 3.9S 

171. 8M 

A1  Q 

Albuquerque,  New  Mexico 

34. 9 N 

106. 5 W 

cue 

Chiang  Mai,  Thailand 

18. 8N 

99 . OF 

CTA 

Charters  Towers,  Australia 

20.  IS 

146. 3M 

l)AV 

Davas,  Philippines 

7.  IN 

125.6M 

ElL 

Eilat,  Israel 

29. 3N 

34. 5M 

FBK 

Fairbanks,  Alaska 

64. 9N 

148. 0W 

KIP 

Kipapa,  Hawaii 

21.4N 

158. 0W 

KON 

Kongberg,  Norway 

59 . 6N 

9.6E 

MAT 

Matsushiro,  Japan 

36. 5N 

138. 3M 

MFP 

Fernando  Poo,  Cameroon 

3.3N 

8.7E 

NT)  I 

New  Delhi,  India 

28. 6N 

77. 2E 

ogd 

Ogdenburg,  New  Jersey 

41. IN 

74. 6W 

PKK 

Pretoria,  South  Africa 

25. 7S 

28. 2E 

I’SC 

Faster  Island,  Pacific  Ocean 

27. IS 

109. 4W 

SBA 

Scott  Base,  Antarctica 

77.9S 

166. 8M 

SJC 

San  Juan,  Puerto  Rico 

18.  IN 

66. 1W 

SOM 

Sombrero,  Chile 

52. 8S 

69 . 2W 

TOL 

Toledo,  Spain 

39. 9 N 

4.0W 

ZLP 

La  Paz,  Bolivia 

16. 3S 

68. 1W 

5-stat ion 

g 1 ob  a 1 

network:  CTA,  MIL,  FBK,  PRM, 

ZLP 

5-stat ion 

northern  hemisphere  network:  CUC,  M 

IF,  FBK,  KIP,  SJC. 

10-stat ion 

global 

network:  ALQ,  CMC.,  CIA,  MIL, 

FBK,  MAT,  PRM,  PSC 

, SBA,  ZLP 

10-stat ion 

northern  hemisphere  network:  ALQ,  Clin , DAV,  Ell,  FBK, 

KIP,  KON, 

MAT,  (XII) , SJC 


20-station  global  network:  all 
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»n.»ski*J  > 


detection  probability  at  a given  ugnitude  and  multiplying  the  number  of 

reported  masked  events  bv  this  faeu.r  to  get  an  estimate  of  the  number  of 
truly  masked  events.) 

11,0  i>Vt,U  Saaiplt>  in  th0  simulation  runs  was  1000,  whether  or  not 

aftershocks  were  included.  Although  the  event  sets  were  completely  different 
both  event  simulation  models,  the  ten  network  runs  in  Table  VI  for  each 
model  Were  made  beginning  with  the  same  random  number;  that  is,  all  ten 
networks  "looked"  at  exactly  the  same  sot  of  1000  events  in  each  case.  The 
results  with  "normal"  thresholds  are  for  when  the  background  noise  values 
Were  set  as  previously  stated;  and  the  "low"  threshold  results  are  for  when 
these  noise  values  are  reduced  by  a factor  of  three,  a reasonable  possibility 

if  Sna11  arrayS  rG(,laced  the  single  seismometers  at  each  site.  Cumulative 

detection  results  have  been  computed  using  both  a = 3.0  and  a M =3.5 

cutoff,  Rocull  that  events  with  M = l n . * ^ i . 

s • the  lowest  magnitude  epicenters 

generated  in  the  simulation;  and  Mg  = 3.5  is  approximately  equivalent  to  a 
“b  * 4,0  eartHqUake  Usin«  the  Ma-“b  f°rnula  «ivan  Previously  in  this  report. 

Ihe  percentages  of  mixed  and  masked  events  in  Table  VI  are  somewhat 
less  than  to  be  expected  with  real  detection  due  to  reasons  previously 
mentioned:  1)  deep  events  are  not  simulated,  2)  the  body-wave  portion  of  the 

seismogram  is  not  included  in  the  simulation,  and  3)  masking  in  the  simulation 
run  is  declared  only  if  the  event  was  detectable  under  normal  conditions. 

Also  the  Percentage  of  mixed  and  masked  events  is  almost  always  higher, 
as  expected,  tor  the  case  where  aftershocks  are  included  in  the  model. 

Those  instances  where  they  are  not  is  due  to  the  fact  that  the  event  sets 
are  completely  different  for  the  two  cases,  and  1000  events  is  not  enough 
to  smooth  out  statistical  fluctuations  of  the  sample.  As  noted  before,  we 
choose  an  aftershock  model  which  we  believed  produced  more  events  closely 
spaced  in  time  and  location  than  real  seismic  activity  and  a main  shock  model 
which  is  completely  random  in  time  and  location;  thus  the  results  for  the 
cases  with  aftershocks  in  Table  VI  should  more  or  less  approach  the  real 
situation.  Anyway,  the  relative  amounts  of  mixing  and  masking  in  the  various 

Inflation  runs  should  be  meaningful  even  if  their  absolute  amount  is  some- 
what in  prror. 
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In  regard  to  the  results  In  Table  VI  then,  we  note  the  following: 

1)  The  station  statistics  (actually  an  average  over  all  stations  in 
the  network  under  consideration)  are  stable.  This  is  as  expect ■ I 
since  all  stations  in  all  five  networks  are  roughly  of  equal 

capability. 

2)  Between  5%  and  8 % of  the  events  with  M 3.5  (n^  4.0)  will  r.ot 

be  clearly  detected  at  a typical  I.PE  long-period  site  due  to 
interference.  (These  percentages  should  be  somewhat  higher  for 
detection  of  real  events.)  If  station  thresholds  were  lowered 

by  .5  Mg,  the  corresponding  percentage  would  be  between  20%  and 
242  while  the  number  of  detections  would  be  roughly  50%  greater. 

3)  Although  a seismic  network  detects  more  events  as  its  size  is 
augmented,  the  number  of  mixed  and  masked  events  is  apparently 
decreased;  this  is  due  to  the  fact  that,  with  more  stations, 
there  are  more  chances  for  events  (provided  they  are  not  after- 
shocks) to  be  separated  sufficiently  in  arrival  time  for  detection. 

4)  There  is  no  significant  difference  between  a global  network  and 
a more  limited  network  in  the  northern  hemisphere  in  regard  to 
overall  detection  thresholds  and  numbers  of  mixed  or  masked  events. 

5)  When  stations  thresholds  are  decreased  bv  0. 5 M , the  number  of 

s 

mixed  and  masked  events  reported  by  a network  does  not  increase 

significantly  for  3.5;  and  the  interference  problem  prevents 

all  of  the  simulation  networks  with  low  thresholds  from  achieving 

100%  detection  down  to  M = 3.5.  However,  for  M 3.0,  lower 

s s — 

thresholds  cause  a considerable  increase  in  network  mixed  and 
masked  events  due  to  the  greater  number  of  events  between  3.0 
and  3.5  (64.5%  of  the  randomly  generated  events  had  3.0  M i.j. 

— c 
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CONCLUSIONS 


I.ong-period  interference  is  a significant  problem  which  is  not  accounted 
lor  in  conventional  station  or  network  capability  estimates,  such  as  by 
program  NETWORTH  at  the  SDAC.  .nterference  depends  on  the  magnitude,  shape, 
and  duration  of  the  seismic  signals,  which  vary  considerably  depending  on 
path  and  station.  Nonetheless,  a reasonably  accurate  model  of  signal  Inter- 
ference can  be  programmed  by  parameterlting  average  coda  shapes  and  durations. 
Such  a program  can  be  calibrated  against  real  data  and  then  utilised  to 
predict  interference  effects  on  proposed  networks.  In  this  way  we  were  able 
to  predict  that  neither  augmenting  a network  of  stations  nor  lowering  the 
thresholds  at  its  stations  will  cause  an  Increase  In  the  number  of  mixed  and 
masked  events  above  a certain  N,  even  though  more  are  detected  above  this 
Ms.  lowering  the  station  thresholds  will  produce  more  mixed  and  masked  events 
as  this  Ms  cutoff  is  decreased  though;  on  the  other  hand  augmenting  the 
network  will  not  cause  a similar  effect  at  these  lower  M values.  Clearly, 
when  Ms  thresholds  are  pushed  lower  and  lower,  there  will  be  more  and  more 
mixed  and  masked  events  to  be  sorted  out  since  the  number  of  events  increases 
with  decreasing  magnitude.  A long-period  system  which  is  designed  to  sort 
out  these  lnterf erring  events  and  which  works  well  at  current  threshold  levels 
may  well  be  overwhelmed  when  thresholds  are  lowered  considerably. 
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APPENDIX 


Normalized  plots  of  amplitudes  of  Rayleigh-wave 
signals  and  coda  recorded  at  LPE  sites 
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